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ABSTRACT 

Mid-infrared atomic and ionic line ratios measured in spectra of pre-main 
sequence stars are sensitive indicators of the hardness of the radiation field im- 
pinging on the disk surface. We present a low-resolution Spitzer IRS search for 
[Ar II] at 6.98 /im, [Ne II] at 12.81 fim, and [Ne III] 15.55 fim lines in 56 tran- 
sitional disks. These objects, characterized by reduced near-infrared but strong 
far-infrared excess emission, are ideal targets to set constraints on the stellar ra- 
diation field onto the disk because their spectra are not contaminated by shock 
emission from jets/outfiows or by molecular emission lines. After demonstrating 
that we can detect [Ne II] lines and recover their fluxes from the low-resolution 
spectra, here we report the first detections of [Ar II] lines towards protoplan- 
etary disks. We did not detect [Ne III] emission in any of our sources. Our 
[Ne II] /[Ne III] line flux ratios combined with literature data suggest that a soft- 
EUV or X-ray spectrum produces these gas lines. Furthermore, the [Ar II] /[Ne II] 
line flux ratios point to a soft X-ray and / or soft-EU V stellar spectrum as the ion- 
ization source of the [Ar II] and [Ne II] emitting layer of the disk. If the soft 
X-ray component dominates over the EUV than we would expect larger photoe- 
vaporation rates hence a reduction of the time available to form planets. 

Subject headings: circumstellar matter - infrared: stars - planetary systems: 
protoplanetary disks - stars: pre-main-sequence 
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Introduction 



Young stars are often surrounded by circumstellar disks containing gas and dust. Most 
of the disk material is in the f orm of gas, while the dust com ponent constitutes only a small 



fraction of the disk mass (e.g. 



Pascucci fc Tachibana 



20101 ). These circumstellar disks are 



planet nurseries; they provide the raw material for planets to form. Thus, understanding 
their evolution, especially the evolution of the gas component, is key to the understanding 
of planet formation. 

The challenge in tracing the gaseous disk component is that the emission and excitation 
of gas species are very sensitive to chemical abundances, disk temperature, gas density and 
radiation field. Hence, multiple transitions at different wavelengths are necessary to probe 
a range of disk radii as well as of disk heights. As an example, the cold (T ^ 50 K) outer 



> 10 AU) disk is usually probed via CO sub-mm pure rotational lines (e.g. 



Dutrey et al. 



20071 ) while observations of CO rovibrational lines in the ne ar-infrared can be u sed to trace 



the hot (~1,000K) gas within ~1 AU from the central star (INajita et al 



2007ah . 



Intermediate disk radii where terrestrial and giant planets form have been probed 
only recently via new mid-infrared (mid-IR) gas lines discovered with the sensitive IRS 
spectrograph on board the Spitzer Space Telescope (e.g. 



2008 



Salyk et al. 



Pascucci et al. 



2007 



Carr &: Najita 



20081 ). One of the gas lines with a high detection rate in Spitzer spectra 



of young stars is the [Ne II] line at 12.81yum. While this line likely prob es shocked 



gas in an outflow from youn g star-disk systems with high accretion rates (IGuedel et al 



2010 



van Boekel et al 



20091 ). the same line is found to trace the hot (>1,000K) disk 



at mosphe re at several AU from sun- like stars t 



outflows (IPascucci fc Sterzik 



2009 



latdo not have fast (i.e . ~100 km/s 



Najita et al. 



2009 



Sacco et al. 



2OI2I). Hence, the 



[Ne II] line at 12.81 fim provides a unique diagnostic to study th e region of the disk th at 



is exposed to high energy stellar photons from the central star (IGlassgold et al. 



2007 
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Gorti fc HoUenbachlboOSl : iHollenbach fc GortilboOQl : lErcolano fc Owenlboioh . Because this 
hot atmosphere can become unbound at large disk radii, [Ne II] hues can be also used 
to test disk dispersal mechanisms such as photoevaporation dr iven by the central star 



( lAlexander 



2008 



Pascucci fc Sterzik 



2009 



Pascucci et al. 



20111 ) ■ 



An interesting class of protoplanetary disks that might reveal how disks d isperse their 



Strom et al. 



19891). These 



gas and dust is that of the so-called transitional disks (e.g. 
objects display small near-IR but large far-IR excesses pointing to an optically thin inner 
cavity within the dust disk, believed to mark the disappearance of the primordial massive 
disk. Such inner cavities could be due to grain growth accompanying terrestrial planet 
formation, photoeva poration, dynam ical clearing by a giant planet or a combination of 



these processes (e.g.. 



Armitage 



20111 ). 



Recently, Spitzer high-resolution mid-infrared spectra of transitional disks have been 
discussed in the literature. In contrast to classical protoplanetary disks whose mid-IR 



spect ra are dominated by molecular emission lines, especially of water (jCarr fc Najita 



20081 ). the spectra of transitional disks pre sent weak molecular lines but strong atomic 



and ionic emission lines (e.g. from [Ne II], 



Najita et al. 



2010 



Pontoppidan et al. 



20101 ). 



In addition, transitional disks do not have signs of outflow activity pointing to a disk 
origin for most of the mid-infrared lines detected in their spectra. In summary, the lack 
of strong molecular lines and outflow activity make transitional disks ideal targets for 
detecting and analyzing additional atomic/ionic lines that could trace disk gas even using 
the low-resolution (hereafter LR with A/AA ~ 60-120) IRS spectrograph on the Spitzer 
Space Telescope. 

A very interesting ionic line that cannot be accessed with the high-resolution (hereafter 



HR with A/AA ~ 600) 



theoretical disk models (IHollenbach Sz Gorti 



RS module is the [Ar II] line at 6.98/im. According to recent 



20091 ). the [Ar II] line luminosity should be 
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similar to that of the [Ne II] line in a disk irradiated by stellar EUV (13.6 eV < hi^ < 100 
eV) and X-ray (0.1-2 keV) photons. In addition, the [Ar II] /[Ne II] and [Ne II] /[Ne III] 
line ratios could help discriminate between cases in which a hard EUV, or a soft EUV or 
X-ray stellar spectrunj^ provide the main source of ionization for the disk atmosphere. The 
[Ar III] line at 8.99 fim could be another important diagnostic of the stellar radiation field 
impinging on the disk and it is often detected in ionized gas surrounding massive stars 



Mart in- Hernandez et al 



2OO2I) 



(H II regions) together with [Ne II] and [Ar II] lines (e.g. 
However, this line falls within the strong, broad, and common 10 fim silicate emission 
featu re also arising from the di sk surface and it is expected to be weaker (by about a factor 



of 5, 



HoUenbach fc Gorti 



20091 ) than the [Ne II] and [Ar II] lines. Therefore, we have not 
attempted to recover this line in these low-resolution spectra. Constraining the main source 
of disk heating and ionization is necessary for properly estimating the rate at which disk gas 
is photoevaporated and hence the predicted disk lifetime. For instance, if photoevaporation 
is induced solely by EUV photoi is then photoevaporat ion rates are expected to be very low. 



of the order of ~ 10 ^^MqJt ^ ( lAlexander et al. 



20041 ). Stellar X-rays can penetrate much 



larger columns of gas than EUV photons resulting in higher photoevap oration rates, up to 



two orders of ma gnitude higher than EUV-induced photoevaporation ( iGorti et al. 



2009 



Owen et al. 



20101 ). In this paper, we conduct a Spitzer-hased search for [Ar II] 6.98 /im. 



[Ne II] 12.81 /im, and [Ne HI] 15.5 fim lines in a large sample of transitional disks. The 
resulting line ratios are compared to models of disks irradiated by high-energy stellar EUV 



^A "hard" spectrum is that in which there are more photons at higher energies beyond a 
certain energy value. A hard X-ray spectrum is that in which there is substantial contribution 
from photons at and beyond 1 keV, while a hard EUV spectrum is that in which there are 
more photons with energies above 41 eV than below this value (for instance the often assumed 
power law for the EUV luminosity but also a black body with temperature >35,000K). 
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and X-ray photons to constrain the radiation field reaching the disk atmosphere. 



2. Sample Selection and Data Reduction 

We collected a comprehensive sample of disks from the literature that have lower near-, 
to mid-IR excess emission than classical protoplanetary disks. This sample includes classical 
transitional disks, cold disks, pre-transitional disks, and anemic diskqj. In the following 
we will use the term transitional disks to refer to this sample of objects. We collected all 
available archival Spitzer IRS spectra for our sample. Since our main goal was to detect the 
[Ar II] line at 6.98 fim and the [Ne II] line at 12.81 /im, we focused on the Short Low (SL) 
module of the Spitzer IRS spectrograph, which covers the wavelength region between 5.2 
and 14.5 /xm. However, when available we also reduced Long Low (LL) spectra to search 
for the [Ne III] line at 15.55 fim. Based on the post-BCD processing, we excluded spectra 
from the sample that had poor signal-to-noise, or were contaminated by the spectrum of a 
nearby source, or the source was not detected in the SL module. Our final sample contains 
56 objects covering an age range of 1-10 Myr, which are listed in Table [1] together with 
their classification and observation log. 

Our data reduction started from the droopres intermediate data products obtained 
via processing of the raw IRS LR data with the SSC pipeline S18.7.0. Further processing 
was done with the spectral reduction tools developed for the Formation and Evolution o f 



Plane tary Systems (FEPS) Spitzer Legacy program (IMeyer et al. 



2004 



Bouwman et al. 



20081 ). In brief, the pairs of imaged spectra acquired along the spatial direction of the 
slit were subtracted from each other in order to correct for the background emission and 
stray light. Then, bad pixels were replaced with an interpolated value obtained from 



^http:/ /www. ipac.caltech.edu/spitzer2008/talks/Diskionary.pdf 
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neighboring good pixels. Finally, the spectra were extracted from the background-subtracted 
pixel-corrected images using a 6.0 pixel and 5.0 pixel fixed-width aperture in the spatial 
dimension for the observations with the SL and LL modules, respectively. The spectra 
were calibrated using a spectral response function derived from multiple IRS spectra of 
the calibration star t]^ Doradus and a MARCS stellar model provided by the Spitzer 
Science Center. Internal uncertainties per pixel were estimated by the FEPS pipeline as 
the standard deviati on of the repeated spe ctral observations (i.e. for all slit positions and 



cycles, see details in 



Carpenter et al. 



20081) 



In order to identify [Ne II] , [Ne III] and [Ar II] lines and to measure line fiuxes, we used 
a robust, Monte-Carlo based Gaussian fitting method with IDL. In the case of [Ne II] lines, 
we fitted simultaneously a Gaussian to the line and a local parabola to the continuum (see 
Figured]). We allowed at most a centroid shift of 0.013 fim (roughly 20% of the wavelength 
spacing) to account for the wavelength calibration uncertaintie in the SLl module (IRS 
Instrument Handbook v. 4 If]). Because at the low-resolution of Spitzer (>2,000km/s) lines 
from disks and even from jets/outfiows are unresolved we further fixed the Gaussian width 
to 0.043 ± 0.008 /im, which is the value and uncertainty we determined from fitting 7 
strong and unresolved [Ne II] lines in different supernova remnant LR Spitzer spectra. 
For the Monte-Carlo approach, we generated 1000 spectra within the range of a Gaussian 
distribution with a sigma equal to the uncertainty at each wavelength. For all of the spectra 
we repeated the fitting method and computed the area of the Gaussian. We used the mean 



of the line fiuxes as final line fiux and the standard deviation o 
the uncertainty over the final line fiux (see 



Pascucci et al. 



2008 



and 



generated line 


iuxes 


Banzatti et al. 


2012 



for 



a similar method applied to silicate emission features and gas lines). 

In the case of the [Ar II] line the adjacent intervals were too noisy for a simple Gaussian 



^http:/ /irsa.ipac. caltech.edu/data/SPITZER/docs/irs/irsinstrumenthandbook/ 
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fit. Thus, here we first subtracted a continuum with the following approach. We produced 
a median spectrum for each object which had multiple IRS LR spectra, in order to identify 
common lines and bands in the SL2 module (5.13-7.60 /im). Then, excluding these lines 
and bands we fitted a robust smooth curve to this module. We added in quadrature the 
continuum subtraction uncertainty to the original uncertainty on each datapoints. After, 
we followed the method described above for the [Ne II] lines to determine the line flux and 
uncertainty (see Figure |2]). Here, the allowed centroid shifts were 0.008 /im (roughly 20% of 
the wavelength spacing) which is the wavelength calibration uncertainty for the SL2 module 
(IRS Instrument Handbook v.4). The Gaussian width parameter was fixed to 0.0241 ± 
0.0041 fim since this is the width and uncertainty of unresolved [Ar II] lines we measured 
from supernova remnant spectra as done for the [Ne II] lines. 

When no line was detected, we computed 3a upper limits with the following equation: 

F„p = SRdxVN (1) 

where R is the RMS of the data points around the line, and the line flux is computed over 
the resolution element {N = 2 pixels), assuming the noise is uncorrelated. Here, dx is 
the wavelength of the given line divided by the resolution (fixed to 120). Along with the 
measured fluxes, the upper limits are in Table El 

As a check to the absolute flux calibration, we integrated the IRS spectrum in the 
IRAC 8 fim filter band and compared the integrated IRS flux to non-contemporaneous 
IRAC fluxes from the literature. Among the 31 objects that have IRAC fluxes, 42% showed 
less than 10% difference and 55% showed less than 20% difference. The absolute flux 
calibration accura cy is ^ 10%, beca. use we did not correct the flux for the location of the 



source in the slit (jSwain et al. 



20081). That about half of the objects p resent difference s 



i n flux es larger than 20% is likely due to intrinsic infrared variability. 



Espaillat et al 



( I2OIII ) found that the infrared spectral energy distributions (SEDs) of 8 out of 14 pre- 
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and classical-transitional disks vary by more than 20% on a timescale of years, consistent 
with our IRS-IRAC flux comparison. Because the [Ar II] and [Ne II] lines fall on different 
segments of the SL module, we performed additional tests to make sure that mispointing 
would not affect the line flux ratio. First of all, we checked that the SL and LL (when 
available) spectra match to within 10%. In addition, based on the difference in fluxes at 
overlapping wavelengths between the SLl, SL2, and SL3 modules we found that we would 
need to apply only small offsets in the source location within the slit to obtain the best 
match between these modules. These offsets do not affect line fluxes within the error bars 
we report here. As a demonstration of this, we applied these small offsets to four objects 



that have both [Ar II] and [Js 
[PZ99]J160421. 7-213028, see 



e II] emissions in their spe ctra (TW Hya, CS Cha, 16201-2410, 



Sicilia-Aguilar et al. 



2009 



for the procedure). We could verify 



that the line fluxes (and hence their ratios) reported without applying any offset is the same 
as with applying the offsets within the uncertainties that we assign. 



3. Results 

Although the LR module has a spectral resolution that is 5 times lower than that 



disks as demonstrated by 



Pascucci et al. 


( 


2009 


) 


2; study by 


Teske et al. 



or the HCN and C2H2 rotation-vibration 



( 1201 ll ). Here, we demonstrate that even 



atomic lines can be recovered in LR spe ctra of transitiona l disks. Figure [3] presents the 



TWHya high-resolution SH2 spectrum ( iNajita et al. 



20101 ) and our LR spectrum around 



the [Ne II] 12.81 fim emission line of the same object. The comparison of the two spectra 
shows that the strongest lines in the HR module, the H I (7-6) line at 12.37 /im and the 
[Ne II] line at 12.81/im, are also detected in the SL module. In addition, we show that 
the flux can be also recovered in LR spectra. The [Ne II] line flux measured in the LR 
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spectrum is (5.8 ± 0.8) x 10 erg s ^ cm ^, which agrees within the estima ted error bar 



with the flux from the HR SH2 spectrum: (5.56±Qg^) x 10 erg s ^ cm ^ (INajita et al 



2010l ). In addition to TW Hya, we also examined 9 other objects from our sample, which 



have availab 


e [Ne 


Lahuis et al. 


2007 



I] line fluxes in the literature 



Baldovin-Saavedra et al 



( iGuedel et al. 



2010 



Pascucci et al. 



2007 



201ll ). We compared the fluxes from HR 



spectra to our LR flux values and we found that the fluxes are in good agreement within 
the associated error bars (see Figure H]). 

We report here the flrst detections of [Ar II] lines in protoplanetary disks (Fig. 
For the complete transitional disk sample, signiflcance of each line detection is computed 
by dividing the line flux by its formal uncertainty (see Table [2]). Two objects feature this 
line at a level >3 a and several others slightly below this threshold. Our detections prove 
that, in spite of their lower spectral resolution, the Spitzer LR spectra can be used to 



detect \At III 


ines. 


Ne III lines are o" 


ten detected in Spitzer HR spectra of disks 


e.g. 


Pascucci et al. 


2007; 


Lahuis et al. 


2007; 


Guedel et al. 


2010; 


Baldovin-Saavedra et al. 


2011 



see Table [3]). We show here that detections are common even in LR spectra of transitional 
disks. In our sample [Ne II] emission appears in 17 objects with ^ 3a level (see Fig. [1]). 
Altogether one object from the total of 56 has both [Ar II] and [Ne II] gas lines in the same 
spectrum, counting only the ^ 3a detections. In addition, one object shows only the [Ar 11] 
line while 16 have only the [Ne II] line. We did not detect the [Ne HI] line at 15.5 fim in 
any of our objects, in agreement with the very few [Ne HI] line detections in HR spectra 
in the literature. This shows that the [Ne HI] line is typically weaker than the [Ne II] line 
even in transitional disks. 

Because multiple Spitzer/IRS data reduct ions have been perforrn ed - e.g. Cornell 



Atlas of Spitzer/IRS Sources (hereafter Cassis, iLebouteiller et al. 



20 111 ) and SSC Enhanced 



Products 







When multiple observations of the same source were available, we found variability in 
the continuum and also in the intensity of the gas lines. There were multiple observations 
when we did not find a line that was however detected in most of the other spectra 
(e.g. GM Aur, see Fig. [T]). We attribute variations by more than 20% to intrinsic 



source variability which was alread y observed in the mid- 



featu res of circumstellar disks (e.g. 



MuzeroUe et al 



2009 



R continuum and dust emission 



Bary et al 



2009 



Espaillat et al. 



201ll ). Infrared va riability is more like ly a common characteristic of disks, especially of 



transitional disks ( lEspaillat et al. 



201l[ ). We emphasize however, that this variability does 



not significantly affect our analysis which focuses on flux ratios of lines belonging to the 
same mid-infrared spectrum. 

In the simplest scenario in which ionization is carried out by EUV photons, the 
line flux ratio from two successive stages of ionization of a given element can be used 
to measure the slope of the sp ectral energy distribution (SED) of the ionizing field (e.g. 



Martin-Hernandez et al 



2OO2I ). X-ray irradiation and ionization also produce different line 



ratios depending on the hardness of the X-ray spectrum. In the following we compute the 
[Ne II] /[Ne III] and [Ne II] /[Ar II] line flux ratios and analyze their distribution. The results 
will be compared with models from the literature in order to gain insight into the radiation 
field impinging on the disk surface. 



%ttp:/ /irsa.ipac. caltech.edu/data/SPITZER/docs/irs/irsinstrumenthandbook/84/ - we 
compared our reduced spectra with the products from those different pipelines. We down- 
loaded the same observations (when they were available) and run our line detection routines 
on these spectra. We found, that the detections were the same (for [Ne II] and [Ar II] lines 
as well) regardless of the used reduction pipeline. Usually we got similar SNRs, see examples 
on Fig. [2]). 
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3.1. [Ne II] /[Ne III] Line Flux Ratios 



Because of the [Ne III] line non-detections in our LR spectra, we can only derive 
lower limits for [Ne II] /[Ne III] line ratios (see left panel on Figure [5]). On Fig. |5] 
we also over-plot model predictions from: a hard EUV spectrum {L^ oc in light 
gray, and a soft EUV (black body e mission at Tcs =30,000 K)/X-ray irradiation in 



dark grav (IHoUenbach &: Gorti 



HoUenbach fc Gortil fl2009h 



20091). Expanding o n the X-ray model predictions by 



Ercolano fc Owen! (120101 ) simulated 3 different X-ray luminosity 



cases for primordial disks and various hole sizes for transitional disks. We note that in these 
models the hole is both in the gas and in the dust distribution while most transitional disks 
are not completely empty of gas in the dust hole as evinced by their on-going accretion. 
These model predictions are also plotted on Fig. [5] with light green slanted stripes for 
various X-ray luminosity models and light blue vertical stripes for various hole size models. 
In summary, all model predictions point to a [Ne II] /[Ne III] line ratio > 1 for an impinging 
soft EUV spectrum and/or X-ray (hard or soft) spectrum. 

Since all our data points are lower limits, they cannot unequivocally discriminate 
between hard EUV and X-ray/soft EUV models. However, several of our lower limits lie 
above the hard EUV limit, suggesting the dominance of the X-ray/soft EUV spectrum. 
In order to further analyze this question we collected [Ne II] /[Ne III] line flux ratios from 
the literature computed from HR IRS spectra (Table [3]). Even with the HR module most 
[Ne III] lines are not det ected. The except ions are the [Ne III] lines from Sz 102 , an 



object with a known jet (iLahuis et al. 



20071 ). WL5/GY246 flFlaccomio et al. 



possible detection of this line from TW Hya ( INajita et al. 



20091), and a 



2010l ). These lower limits and 



measurements from the literature on the right panel of Fig. [5] clearly favor an X-ray/soft 
EUV ionizing spectrum impinging on the disk. All [Ne II] /[Ne III] lower limits from the HR 
spectra are > 1 and thus exclude that a hard EUV spectrum ionize neon atoms at the disk 
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surf ace. This result is in agreement with pr evious HR studies of smaller samples of disks 



(e.g. 



Pascucci et al. 



2007 



Lahuis et al 



20071 ) 



3.2. [Ne II] /[Ar II] Line Flux Ratios 

Figure |6] shows the distribution of the computed [Ne II] /[Ar II] line flux ratios for our 
sample of transitional disks. The dataset contains measured values, lower limits and also 



one upper limit. On t he same figure we a. 



model predictions by 



HoUenbach fc Gorti 



s o ove r-plot the hard X-ray and soft X-rays/EUV 



( I2OO9I ) with darker gray and vertical light gray 



stripes, respectively. Their model uses solar Ne/Ar elemental abundance ratio. Due to 
the dominance of lower limits we cannot discriminate between the models. However the 
distribution of the data, in particular the measured flux ratios, cluster closer to the ~ 1 
which points to the soft X-ray/EUV case. 

We also computed the median of continuum-subtracted and distance-corrected spectra 
for [Ne II] line detections and found a median flux value of 1.49 xlO~^^ erg s~^ cm~^ at the 
distance of the Taurus star-forming region (140 pc). In the case of the [Ar II] line, where we 
have only 2 firm detections, the fluxes at the Taurus distance are 1.81 and 5.90 xlO~^^ erg 

— 1 —2 

s cm . 

Our sample is dominated by [Ne II] /[Ar II] line ratio lower limits, hence we perform 
further analysis to identify what is the typical line flux ratio for the sample. In order to 



compute the representative line ratio values for the distribution, we app^ 
with the software package "Astronomy SURVival A nalysis" (ASURV 



1990 



ied survival analy sis 



Lavalley et al. 



1992 



Feigelson fc Nelson 



Isobe &: Feigelson 



19851 ). We chose this statistical method due 



to the many lower limits in the flux ratio data. Survival analysis can provide a cumulative 
distribution function of a given sample containing measured data points and censored 
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data points as well. From the distribution, this method can also estimate a mean value 
of the distribution. Among the ASURV different subroutines, we applied a Kaplan-Meier 
estimator, which is a self-consistent, generalized maximum-likelihood estimator for the 
population from which the sample was drawn. The ASURV code can handle only right 
or left censored data (only lower or upper limits respectively), but cannot do both at the 
same time. Therefore we applied the analysis on [Ne II] /[Ar II] line ratio lower limits and 
measured data points. As a result, the mean is equal to 1.20 ± 0.06. This value is below 
the hard X-ray stellar spectrum model hinting to the soft X-ray and EUV spectrum model 
Although, given that the majority of the [Ne II] /[Ar II] line ratios are lower limits, one 
cannot unequivocally rule out the hard X-ray model either. 



4. Discussion 



Emission from ionized atoms, such as Ne"*", Ne"*'"*', and Ar"*", can be produced 
either in the interstellar medium gas that is shocked by protostellar winds or in the 



l ow-mass-disk surface layer 



( iHollenbach Sz Gorti 



2009 



by gas that is heated and ionized by ste 



Ercolano fc Owen 



2010 



Glassgold et al 



ar hi gh energy photons 



20071). Aspect acular 



van Boekel et al. 



(120091 ) spatially 



example of the first case is the T Tauri triple system where 
and spectrally resolved the strong [Ne II] emission detected in Spitzer spectra. However, 
transitional disks, such as those in our sample, have no obvious signs of outflow activity 
suggesting that the unresolved emission lines in Spitzer spectra trace the disk surface. This 
has been confirmed at least in seven transition al disks observed at high spectral and spati al 



2009 



Sacco et al. 



20121 ). 



resolution with the VLT/VISIR spectrograph (jPascucci fc SterzikI 
The [Ne II] FWHMs are found to be relatively narrow (15-40 km/s) and to peak close to 
(but not exactly at) the stellar velocity. An outflow origin for this line would result in much 



broader profiles (~100km/s) as well as more blueshifted peaks 



-100 km/s). Further 
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evidence th at the [Ne II] e missio n in transitional disks arises primarily from the disk is 



provided by 



Guedel et al. 



( l2010l ). They find that known outfiow sources are remarkably 
separated from the other sources by having [Ne II] luminosities 1 to 2 orders of magnitude 
higher. These findings strongly support a disk origin for the [Ne II] emission detected 
toward transitional objects, and very likely for the [Ar II] lines as well. In the specific 
case of the transitional disk around TW Hya, the spectrally resolved [Ne II] line profiles 



demonstrate that most of t 



radius fjPascucci et al. 



le [Ne II] emission arises from the disk beyond the dust inner 



20111). 



We turn now to model predictions of disks irradiated by EUV and X-ray photons to 
understa nd what our observa t ions a re telling us about the radiation impinging on the disk 



surface. iHollenbach &: Gortil (j2009l ) have shown analytically and with numerical models 
that a hard EUV (J-ieuv{u) oc spectrum produces more [Ne III] than [Ne II] emission. 
This stems from the fact that there are more high-energy (around 40 eV) than low-energy 
(around 21 eV) photons ionizing Ne atoms and Ne"*" ions, as also found in HII regions 
when the stellar temperature of the massive star is >35,000K. A higher [Ne III]/[Ne II] 
line fiux ratio is clearly in contrast with observed ratios (see Sect. 13.11) thus excluding a 
hard EUV spectrum as the main source of ionization for neon atoms in the atmosphere 
of transitional disks. Because the detected [Ne II] lines are at least ~10 times stronger 
than the [Ne III] lines, if they trace the EUV layer then the EUV spectrum impinging on 
the disk should be that created by a black body with effective temperature 30,000 K (the 
spe ctrum drops sharply frorn 21 eV to 40 eV which results in more low-energy photons. 



see 



HoUenbach fc Gorti 



20091 ). A disk irradiated so lely by X-rays also produces more 



[Ne II] than [Ne III] emission ( iGlassgold et al. 
reactions of Ne"^"*^ with atomic hydrogen. This [Ne 
agreement with observations. Further modeling by 



20071 ). because of the rapid charge exchange 
I] over [Ne III] do i ninan ce is also in 



Ercolano fc OwenI (120101 ) including 



both X-rays and EUV shows that this ratio is sensitive to the X-ray luminosity as well as 
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to the disk structure, transitional versus classical diskqj. According to their model, the 
[Ne II] /[Ne III] line flux ratio ranges from ~ 3 to ~ 5 in disks with inner holes of different 
sizes from 8.5 AU to 30.5 AU. Moreover, this ratio ranges from ~ 5 to ~ 9 in primordial 
disks depending on various X-ray luminosities. In summary, the observed [Ne II] /[Ne III] 
line ratio could be produced either by a soft EUV, a soft X-ray, or a hard X-ray spectrum. 

The first ionization potential of argon is 15.76 eV, lower than that of neon. Thus, 
the intensity of [Ar II] lines allows one to constrain the SED of the ion izing field at lower 



HoUenbach fc Gorti 



energi es than what is possible with the [Ne II] and [Ne III] transitions. 
(120091 ) predict the [Ar II] 6.98 /im line to be one of the strongest forbidden lines arising 
from the hot disk atmosphere. More specifically, they predict a strength very similar to 
that of the [Ne II] 12.81 fim line (assuming solar Ne/Ar elemental abundance ratio) when 
the heating is by EUV or soft X-rays and a factor of ~2.5 weaker line when the heating 
is by hard X-rays (see also their Fig 4). This is mainly due to the fact that the X-ray 
heated gas is at a lower temperature, closer to that of the [Ne II] transition than to the 
much higher [Ar II] line. The line flux ratios measured in this work point to a soft EUV 
and/or soft X-ray spectrum producing the [Ne II] and [Ar II] lines in transitional disks. 
These soft emission spectra, in the EUV as well as in the X-ray regimes, are thought to 
be connected to the accretion of disk gas onto the central star. For instance. X-rays could 



arise in accretion s 
(e.g. 



Kastner et al. 



like spectrum (IRibas et al. 



locks along the funnel connecting the circumstellar disk to the star 



20021). Chro mospheric EUV emission is known to have a power-law 



20051 ) ■ thus it could contribute only marginally to the observed 
[Ne II] and [Ar II] emission lines. Based on the [Ne II] li ne median flux value and the 



[Ar II] line flux range reported in Sect. 13.21 and eq. 8 from iHollenbach fc Gortil ( l2009l ) , we 



^A disk that is irradiated by EUV and X-rays will have a fully ionized Hll-like region at 
10,000 K on top of a hot, more neutral X-ray (~1 000 K) layer. 
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also calculate that ~10^^ photons/s are necessary to reproduce these fluxes for a EUV-only 
irradiation. If X-rays contribute or dominate, this number should be taken as an upper 
limit. We note that independent calculations of the stellar ionizing flux produced by 
young stars range from ~10'^°-10^^ photo ns /s for the transitional disk around TW Hya 



flHerczeg et al 



2007 



Pascucci et al. 



stars (Alexander et al. 



20121) and to 10^^-10^^ photons/s for classical T Tauri 
20051 ). Hence, our results are in line with an [Ar II] emission also 
produced by the disk surface. The X-ray component that might dominate disk irradiation 
is especially exciting because X-rays can heat and ionize larger colu i nns of gas than EUV 



and 



]hus drive larger photoevaporative winds (lErcolano et al 



2009 



Gorti fc Hollenbach 



20091 ). Models that include X-ray irradiation of the disk surface predict up to t wo orders o: 



magn itude higher mass loss rates in comparison to EUV-only irradiated disks (jOwen et al. 
2010l ). Such high photoevaporation rates could dominate the evolution and dispersal 
of protoplanetary disks. Ongoing photoevaporation has been detected already in a few 



transitional disks (jPascucci Sz Sterzik 



2009 



Pascucci et al 



2011 



Sacco et al 



2OI2I) 



More [Ar II] line detections are clearly necessary to expand upon these flrst and 
exploratory results. This study will be important for planning future and more sensitive 
MIR observations of disks. The spectrometer MIRI on the James Webb Space Telescope 
will be able to easily detect the ionic lines discussed here. According to current estimates 
MIRI will be able to detect a line flux of 1x10"^^ W/m^ at 7 /im with a signal-to-noise of 
10 in just 100 sec. This flux is about 60 times lower than our [Ar II] line fluxes corrected 
to the 140 pc distance of the Taurus star-forming region. With more [Ar II] detections it 
will be possible not only to better measure the [Ne II] /[Ar II] flux ratio in transitional disks 
but also to explore how line fluxes evolve with stellar and disk evolutionary stage thus 
constraining the time evolution of photoevaporation rates. Conflrming that [Ar II] lines 
really trace photoevaporating gas requires much higher spectral and spatial resolution than 
that available with MIRI (only ~100km/s). EXES on SOFIA is perfectly suited for this 
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task. Although much less sensitive than MIRI, it is still sensitive enough to detect fluxes 
similar to our [Ar II] line fluxes, at the distance of Taurus. Its high resolution mode will 
provide up to R~120,000 or 2.5 km/s around 7 fim enabling to spectrally resolve [Ar II] 
lines if they trace photoevaporating gas at tens of AU from the central star as predicted. 
The 3 signal-to-noise level will be reachable in 140 minutes for a resolved (10 bins) line. Our 
transitional disks with [Ar II] detections are the best targets for high-resolution follow-up 
studies of the gas kinematics with SOFIA. 



5. Summary 

In this paper we present [Ar II] and [Ne II] line detections in low-resolution (LR) 
Spitzer/mS spectra. From the [Ne II] /[Ar II] and [Ne II] /[Ne III] line flux ratios, we 
investigate whether the ionization of the [Ne II], [Ne III] and [Ar II] emitting layer is 
mainly due to soft/hard X-rays or EUV photons impinging on the disk. Our results can be 
summarized as follows: 

1. We report the flrst detections of [Ar II] lines in protoplanetary disks. We detected this 
line in 2 sources at a level >3 a, altogether 4% of the sample show [Ar II] detections. 

2. We also detected [Ne II] lines for the flrst time in LR IRS spectra. Our 17 detections 
with at least 3a level account for 30% of the objects in our sample. 

3. Our [Ne II]/ [Ne III] line ratio when combined with literature data excludes that 
the layer emitting [Ne II] and [Ne III] is mostly ionized by a hard stellar EUV 
(LiEuviu) oc u^^) spectrum. 

4. The [Ar II] /[Ne II] line flux ratios are dominated by lower limits, thus one cannot 
distinguish unambiguously between the two models. However, the distribution of the 
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upper/lower limits, as well as the measured values, seem to favor the soft X-rays/EUV 
stellar spectrum rather than a hard X-ray spectrum reaching the [Ar II] and [Ne II] 
emitting layer of the disk. Clearly, more [Ar II] line detections are needed to better 
constrain photoevaporative disk models. 



A dominance of the soft X-ray component would point to larger photoevaporation rates 



than when photoevaporation is solely driven by 



and extension of dust gaps in transitional disks flOwen et al 



UV photons, wh ich influences the lifetime 



201l|). 
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Fig. 1. — The [Ne II] line detections 3a) where green curves indicate the fitted Gaussians 
and orange curves mark the continuum. The corresponding SNR is plotted on the top left 
corner. 
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Fig. 2. — The two [Ar II] line detections with < 3cr level. For comparison, we plotted the 
same spectra reduced with different pipelines: SSC Enhanced data product in the middle, 
Cassis spectra at the right (see text). The green curves indicate the fitted Gaussians to the 
continuum-subtracted residuals. The corresponding SNR is plotted on the top left corner. 
The line detections are always present in the spectra regardless the used reduction pipelines. 
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Fig. 3. — Comparison of the high-resolution ( iNajita et al. 
spectrum (this paper) of TW Hya. The strongest lines are detected even in the low-resolution 
module. We note that the two spectra were obtained at different epochs. Hence the different 



line ra tio of the HI 
out in lNajita et al. 



7-6) a nd Ne II lines could be due to source variability as already pointed 
( I2OIOI ) for this source based on a second epoch high-resolution spectrum. 
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Fig. 4. — Comparison of the high-resolution (blue, iGuedel et al.l l2010t iNajita et al.l l2010l : 
Baldovin-Saavedra et al.l l201ll ) and the low- resolution [Ne II] fluxes (black). The low- 
resolution line flux values are in good agreement with the high-resolution fluxes. 
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Fig. 5. — Left: [Ne II]/[Ne III] line ratios for each of the objects with [Ne II] emission. We never detected [Ne III] lines 
in our spectra, so all poi nts are lower limits. Right: literature data on [Ne II] /[Ne III] line flux ratios. With darker gray 



background we indicate 



Hollenbach fc Gorti 



(j2009[ ) model for an X-ray/soft EUV dominated stellar spectrum, while 
lighter gray is the hard EUV case. Note that the latter model is less likely given the distribution of lower limits. Light 
green slanted stripes mark Ercolano & Owen (2010) model predictions for various X-ray luminosities, while light blue 
straight stripes are their transitional disk models with various hole sizes. Due to the lower limits, one cannot distinguish 
between these different model predictions. 
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Fig. 6.— [Ne II1/[Ar III 



ine flux ratios for eacli object. We plot model predictions from 
HoUenbach &: Gortil (j2009l ). Light gray vertical stripes show the soft X-ray or EUV model 
predictions, while darker gray stripes indicate the hard X-ray model. Because of the majority 
of datapoints is lower limits, we cannot rule out the hard X-ray model, but the distribution 



hints toward to the soft X-ray or EUV model. 



Table 1. The Transitional Disk Sample 



ID 


Source 


Position {a, 5) 
(J2000) 


Classification 
Criteria 


Ref. 


Exposure Time 
(Ramp X Cycle) 


AORKey 


1 


LAL 31 


03 29 29.25 +31 18 34.7 


a 


1 


14xl3|14x3 


19053312 


2 


SSTc2d 034227.1+314433 


03 42 27.12 +31 44 32.9 


b 


2 


14xl7|30x4 


19053568 


3 


LRLL 67 


03 43 44.61 +32 08 17.7 


a 


1 


60x2|14xl 


22964480 












60x2|14xl 


22970624 












60x2|14xl 


22970112 












60x2|14xl 


22970368 












14x3|6x2 


16267264 


4 


IC348-31 


03 44 18.16 +32 04 57.0 


a* ,c* 


3 


6x2 


16269056 












14x2|6xl 


22963712 












14x2|6xl 


22968576 












14x2|6xl 


22968832 












14x2|6xl 


22969088 


5 


IC348-44 


03 44 22.57 +32 01 53.6 


a* ,c* 


3 


60x2|6x2 


16266752 


6 


LRLL 97 


03 44 25.54 +32 06 17.1 


a 


1 


14x3|30x2 


16267008 


7 


2MASS J03443481+3156552 


03 44 34.81 +31 56 55.2 


b 


2 


14x24|14x34 


19053824 


8 


LRLL 58 


03 44 38.54 +32 08 00.6 


a 


1 


60x2|30x3 


16755456 












14x3|14x2|14xl 


22964224 












14x3|14x2|14xl 


22966016 












60x2 none 


22966272 


9 


IC348-133 


03 44 41.73 +32 12 02.2 


a* ,c* 


3 


60x2|6x2 


16266496 


10 


LRLL 1679 


03 44 52.07 +31 58 25.5 


a 


1 


60x4|30x8 


22971136 












60x3|60x4|30x4 


22970880 












60x4|30x8 


22971392 












60x3|30x4 


22964992 


11 


IC348-114 


03 44 56.14 +32 09 15.2 


a*,c* 


3 


14x2|6xl 


22968832 












6x2 


16269056 












14x2|6xl 


22968576 












14x2|6xl 


22969088 












14x2|6xl 


22963712 


12 


LkHo 330 


03 45 48.28 +32 24 11.8 


b 


4 


14x1 


5634816 


13 


CX Tau 


04 14 47.86 +26 48 11.0 


c 


5 


6X1 


3534592 



Table 1 — Continued 



ID 


Source 


Position {a, 5) 
(J2000) 


Classification 
Criteria 


Ref. 


Exposure Time 
(Ramp X Cycle) 


AORKey 


14 


FO Tau A 


04 14 49.28 +28 12 30.5 


c 


5 


6x1 


3534592 


15 


FQ Tau 


04 19 12.81 +28 29 33.0 


c 


5 


14x2|14xl 


3551744 


16 


UX Tau A 


04 30 03.99 +18 13 49.3 


a 


6 


6x1 


26140928 












6x1 


3536384 












6x1 


27187456 


17 


DM Tau 


04 33 48.71 +18 10 09.9 


a 


6 


60x4|14xl 


16346624 












6x4|none 


15117824 












6x4|none 


15117312 












6x4 none 


15117568 












6x4 none 


19489024 












6x4 none 


19488768 












6x4 none 


19488512 












6x4 none 


15110912 












6x1 


3536384 












6x3|6xl 


26141952 












6x3|6xl 


27184640 


18 


DN Tau 


04 35 27.37 +24 14 58.9 


c 


5 


6x1 


3537152 


19 


LkCa 15 


04 39 17.79 +22 21 03.4 


b 


7 


6x1 


3537664 












6x1 


26140672 












6x1 


27186176 


20 


CoKu Tau/4 


04 41 16.81 +28 40 00.0 


a 


6 


14x2|6xl 


3548416 


21 


GO Tau 


04 43 03.09 +25 20 18.7 


c 


5 


6xl|14xl 


3548928 


22 


GM Aur 


04 55 10.98 +30 21 59.5 


a 


6 


6x4|none 


15118848 












6x3|6xl 


27186688 












6x1 


3538944 












6x4|none 


15119104 












6x4|none 


15119360 












6x4 none 


19483904 












6x4 none 


19483648 












6x4 none 


19483136 












6x4 none 


15111168 



Table 1 — Continued 



ID 


Source 


Position (q, S) 
(J2000) 


Classification 
Criteria 


Ref. 


Exposure Time 
(Ramp X Cycle) 


AORKey 












6x3|6xl 


26141696 


23 


CVSO 224 


05 25 46.75 +01 43 30.3 


b 


8 


60x2|30x5 


16264960 


24 


FM 177 


05 45 41.94 -00 12 05.3 


a 


1 


240x2|60x3|120x2 


12643072 


25 


FM 281 


05 45 53.11 -00 13 24.9 


a 


1 


240x2|60x3|120x2 


12642816 


26 


FM 326 


05 45 56.31 -1-00 07 08.6 


a 


1 


14x2|6x2 


18738944 


27 


FM 515 


05 46 11.86 -1-00 32 25.9 


a 


1 


14x2|14xl|14xl 


12641792 


28 


FM 618 


05 46 22.43 -00 08 52.6 


a 


1 


14x3|14xl|14xl 


12641536 


29 


FM 856 


05 46 44.84 +00 16 59.8 


a 


1 


14x8|14x2 


18748672 


30 


SZ Cha 


10 58 16.77 -77 17 17.0 


a 


6 


6x1 


12696832 












6x1 


27187968 












6x1 


26142464 


31 


TWHya 


11 01 51.91 -34 42 17.0 


b 


9 


6x1 none 


3571456 


32 


CS Cha 


11 02 24.91 -77 33 35.7 


a 


6 


14xl|6xl 


12695808 












6x3|6xl 


26144000 


33 


T21 


11 06 15.41 -77 21 57.0 


c 


10 


6xl|14xl 


12696320 


34 


CHXR 22E 


11 07 13.30 -77 43 49.8 


a 


11 


60x2|30xl6 


18361344 


35 


T25 


11 07 19.15 -76 03 04.8 


a 


6 


14x1 


12695552 












6x2 


26144256 












6x2 


27185152 


36 


T35 


11 08 39.05 -77 16 04.2 


a 


6 


6x1 


27185664 












6x1 


26143488 


37 


C7-1 


11 09 42.60 -77 25 57.8 


a 


11 


60x2|120x2 


12686336 


38 


T54 


11 12 42.68 -77 22 23.0 


a 


6 


14x1 


12695552 


39 


Sz45 (T56) 


11 17 37.01 -77 04 38.1 


c 


10 


6x1 


26142720 












6x1 


27186432 












6xl|14xl 


12696064 


40 


HD 98800 


11 22 05.30 -24 46 39.3 


b,c 


13 


6x1 none 


3571969 












6x1 none 


3571968 


41 


T Cha 


11 57 13.48 -79 21 31.3 


b 


4 


6x1 none 


12679424 


42 


HD 135344 


15 15 48.94 -37 08 55.8 


b 


4 


6x1 none 


3580672 


43 


Sz 84 


15 58 02.52 -37 36 02.7 


b 


2 


14x1130x1 


5644288 
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Table 1 — Continued 



ID 


Source 


Position (o, <5) 


Classification 


Rof. 


Exposure Time 


AORKey 






(J2000) 


Criteria 




(Ramp X Cycle) 




44 


[PZ99]J160421. 7-213028 


16 04 21.65 -21 30 28.4 


c* 


12 


6x8|14x4|120x6|30x4 


19666432 


45 


ScoPMS 31 


16 06 21.96 -19 28 44.5 


a*,c* 


12 


6x4|14x4|6x4 


17777152 


46 


[PGZ2001]J160959.4-180009 


16 09 59.33 -18 00 09.0 


a*,c* 


12 


6x7|6x4|30x6|6x4 


17779968 


47 


SSTc2d J161029.6-392215 


16 10 29.55 -39 22 14.4 


b 


2 


60x4|30xl0 


19051008 


48 


RX J1615.3-3255 


16 15 20.23 -32 55 05.1 


b 


2 


6x2 


15916800 


49 


16126-2235 


16 15 34.56 -22 42 42.1 


a 


6 


6x1 


12675072 


50 


SSTc2d J162245.4-243124 


16 22 45.40 -24 31 23.9 


b 


2 


14x2|30x2 


15920641 


51 


16201-2410 


16 23 09.24 -24 17 04.7 


a 


6 


6X1 


12699392 


52 


SSTc2d J162506.9-235050 


16 25 06.9 -23 50 50 


b 


2 


14x4|30x2 


19059200 


53 


DoAr 21 


16 26 03.01 -24 23 37.9 


a 


1 


14x3|120x3 


23162624 












6xl|14xl 


12698368 


54 


DoAr 28 


16 26 47.41 -23 14 52.1 


a 


6 


14x2|30x2 


12702976 


55 


SR21 


16 27 10.27 -24 19 12.7 


a 


6 


6x1 


12698880 


56 


RX J1852. 3-3700 


18 52 17.30 -37 00 12.0 


c* 


14 


14x4|14x24 


5200640 



Note. — (a) - objects were defined as transitional disks based on their spectral slopes and near IR excess; (b) - SEP modeling; 



(c) 



l 2O10l); (2) 



Merm et al 



(7) 
(12) 



Es£aillat_etal 



Dahm et al. 



(2010D^31 



Cieza et al 



1 20071 '); (Sl lEspaillat et al 



l20o]); (13) 



Furlan et al 



i2007|); (4) 



20081 '); (9) 



l20o])^(14) 



Brown et al. 



Calvet et al 



I 20071) : (5) 



Pascucci et al. 



1 20021) ; (10) 
1 2003) 



Naiita et al. 


I 2007b 


); (6) 


Furlan et al. 


Kim et al. 




fflOsI); (11) 


Luhman et al. 



200£) 



20081) 
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Table 2. Measured Line Fluxes and 3-Sigma Upper Limits 



ID Name [Ar II] flux and uncertainty [Ne II] flux and uncertainty [Ne III] flux and uncertainty 

(lO^i^'erg s^i cm~2) 



1 


LAL_31 


< 0.49 


< 0.22 


< 1.26 


2 


SSTc2d_J034227.1+314433 


< 0.85 


< 0.18 


< 1.14 


3 


LRLL_67.1 


< 0.62 


0.41 ± 0.07 


< 1.97 


3 


LRLL_67 


< 1.40 


< 0.43 


< 2.81 


3 


LRLL_67.2 


< 0.99 


0.37 ± 0.09 


< 2.95 


3 


LRLL_67.3 


< 1.42 


0.25 ± 0.07 


< 1.83 


3 


LRLL_67.4 


< 0.88 


0.30 ± 0.04 


< 1.69 


4 


IC348-31.1 


< 6.28 


< 1.86 


< 2.57 


4 


IC348-31.2 


< 1.91 


< 1.73 


< 2.73 


4 


IC348-31.3 


< 2.07 


< 1.25 


< 5.09 


4 


IC348-31.4 


< 1.65 


< 1.10 


< 4.33 


4 


IC348-31.5 


< 2.29 


< 0.90 


< 3.14 


5 


IC348-44 


< 0.59 


< 0.19 


< 2.57 


6 


LRLL.97 


< 1.23 


< 1.95 


< 1.25 


7 


2MASS.J03443481+3156552 


< 0.75 


< 0.15 


< 0.77 


8 


LRLL_58.1 


< 1.90 


< 0.95 


< 3.71 


8 


LRLL_58.2 


< 3.44 


< 1.76 


< 1.62 


8 


LRLL_58.3 


< 1.69 


< 1.17 




8 


LRLL_58.4 


< 1.46 


< 1.15 


< 1.80 


9 


IC348-133 


< 1.36 


< 0.44 


< 1.57 


10 


LRLL_1679.1 


< 0.69 


< 0.12 


< 0.54 


10 


LRLL_1679.2 


< 0.92 


< 0.08 


< 0.65 


10 


LRLL_1679.3 


< 0.65 


< 0.13 


< 0.53 


10 


LRLL_1679.4 


< 2.71 


< 0.80 


< 0.54 


11 


IC348-114_1 


< 3.40 


< 1.52 


< 2.17 


11 


IC348-114_2 


< 2.49 


< 1.19 


< 2.42 


11 


IC348-114_3 


< 2.14 


< 1.51 


< 3.14 


11 


IC348-114_4 


< 1.82 


< 0.48 


< 6.29 


11 


IC348-114.5 


< 2.81 


< 0.71 


< 3.84 


12 


LkHa_330 


< 5.82 


< 2.13 


< 2.59 


13 


CX.Tau 


< 8.15 


< 1.41 


< 5.79 
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Table 2 — Continued 



ID Name [Ar II] flux and uncertainty [Ne II] flux and uncertainty [Ne III] flux and uncertainty 

(lO^i'^erg cm-2) 



14 


FO.Tau 


< 7.90 


< 2.71 


< 8.85 


15 


FQ.Tau 


< 1.88 


< 0.53 


< 1.74 


16 


UX.Tau.A.l 


< 6.01 


2.09 ± 0.34 


< 4.98 


16 


UX.Tau.A.2 


< 4.72 


< 3.64 


< 4.54 


16 


UX.Tau^.3 


< 8.41 


< 2.33 


< 3.80 


17 


DM_Tau_10 


< 4.26 


0.79 ± 0.22 


< 3.03 


17 


DM_Tau_ll 


< 4.32 


0.94 ± 0.29 


< 2.76 


17 


DM.Tau.l 


< 2.69 


0.80 ± 0.24 




17 


DM_Tau_2 


< 0.51 


0.58 ± 0.07 


< 1.83 


17 


DM_Tau_3 


< 2.98 


< 1.04 




17 


DM_Tau_4 


< 4.26 


< 0.39 




17 


DM_Tau_5 


< 3.32 


0.90 ± 0.30 




17 


DM_Tau_6 


< 1.96 


< 1.12 




17 


DM.Tau.7 


< 3.77 


< 0.54 




17 


DM_Tau_8 


< 3.68 


< 0.78 




17 


DM_Tau_center 


< 4.60 


< 0.91 


< 3.31 


18 


DN.Tau 


< 2.58 


< 2.59 


< 8.46 


19 


LkCa_15.1 


< 3.81 


0.87 ± 0.21 


< 4.31 


19 


LkCa_15.2 


< 5.09 


< 1.31 


< 4.66 


19 


LkCa_15.3 


< 8.59 


< 1.38 


< 2.59 


20 


CoKu.Tau/4 


< 3.54 


< 0.62 


< 4.67 


21 


GO.Tau 


< 8.85 


< 2.09 


< 2.27 


22 


GM_Aur_10 


< 4.98 


1.68 ± 0.39 


< 4.22 


22 


GM_Aur_12 


< 6.19 


1.44 ± 0.37 


< 6.06 


22 


GM^ur.l 


< 4.96 


1.65 ± 0.32 




22 


GM_Aur_2 


< 4.18 


1.12 ± 0.26 




22 


GMJlut.3 


< 3.15 


< 1.16 




22 


GM_Aur_4 


< 3.61 


< 1.36 




22 


GM^ur.5 


< 2.79 


0.94 ± 0.25 




22 


GM^ur.6 


< 2.48 


1.45 ± 0.27 




22 


GM^ur.8 


< 2.07 


1.10 ± 0.26 
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Table 2 — Continued 



ID Name [Ar II] flux and uncertainty [Nc II] flux and uncertainty [Ne III] flux and uncertainty 

(lO^-'-^crg s^-"- cm^-^) 



22 


GM_Aur_9 


< 5.94 


< 1.98 


< 4.57 


23 


CVS0.224 


< 1.03 


< 0.15 


< 0.55 


24 


FM_177 


< 0.78 


< 0.15 


< 0.37 


25 


FM_281 


< 0.63 


< 0.26 


< 0.41 


26 


FM_326 


< 1.42 


< 0.74 


< 1.69 


27 


FM_515 


< 1.76 


< 1.00 


< 1.33 


28 


FM.618 


< 1.92 


0.99 ± 0.08 


< 1.15 


29 


FM_856 


< 1.60 


< 0.39 


< 2.38 


30 


SZ.Cha.l 


< 7.92 


< 2.21 


< 4.88 


30 


SZ.Cha.2 


< 6.15 


< 1.44 


< 3.77 


30 


SZ.Cha.3 


< 4.58 


2.17 ± 0.52 


< 6.86 


31 


TW.Hya 


< 8.35 


5.81 ± 0.81 




32 


CS.Cha.l 


< 2.94 


3.73 ± 0.33 


< 2.37 


32 


CS.cha.2 


< 5.58 


2.83 ± 0.32 


< 2.42 


33 


T21 


< 7.17 


< 2.58 


< 2.13 


34 


CHXR22E 


< 0.69 


< 0.09 


< 0.27 


35 


T25.1 


< 2.09 


< 0.48 


< 1.26 


35 


T25_2 


4.25 ± 1.08 


< 0.60 


< 1.86 


35 


T25.3 


< 4.08 


< 0.66 


< 1.81 


36 


T35 


< 2.68 


1.32 ± 0.29 


< 2.27 


36 


T35.2 


< 8.31 


1.32 ± 0.35 


< 1.93 


37 


C7-1 


< 0.58 


0.10 ± 0.02 


< 0.40 


38 


T54 


< 3.16 


< 1.02 


< 1.15 


39 


Sz.45.1 


< 6.79 


< 1.60 


< 2.67 


39 


Sz.45.2 


< 5.00 


< 1.50 


< 2.03 


39 


Sz.45.3 


< 5.06 


< 1.59 


< 1.16 


40 


HD_98800.2 


< 10.02 


< 4.33 




40 


HD.98800 


< 4.88 


< 5.13 




41 


T.Cha 


< 10.68 


3.09 ± 0.90 




42 


HD_135344 


< 17.81 


< 6.90 




43 


Sz_84 


< 2.25 


< 0.98 


< 1.01 
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Table 2 — Continued 



ID 


Name 


[Ar II] flux and uncertainty 


[Ne II] flux and uncertainty 
(IQ-i^'erg 


[Ne III] flux and uncertainty 


44 


[PZ99]J160421. 7-213028 


1.69 ± 0.49 


1.78 ± 0.14 


< 0.62 


45 


SCOPMS.031 


< 4.19 


< 2.47 


< 2.14 


46 


[PGZ2001]J160959. 4-180009 


< 3.15 


< 1.13 


< 0.92 


47 


SSTc2d_J161029.6-392215 


< 0.30 


< 0.07 


< 0.89 


48 


RXJ1615.3-3255 


< 7.69 


2.97 ± 0.51 


< 2.87 


49 


16126-2235 


< 11.74 


< 1.32 


< 3.65 


50 


SSTc2d_J162245.4- 243124 


< 1.18 


< 0.72 


< 2.04 


51 


16201-2410 


< 10.90 


< 1.53 


< 5.85 


52 


SSTc2d.J162506.9-235050 


< 2.24 


0.73 ± 0.15 


< 1.38 


53 


DoAr_21_l 


< 18.39 


< 23.22 


< 17.66 


53 


DoAr_21_2 


< 17.06 


15.48 ± 3.01 


< 10.39 


54 


DoAr_28 


< 2.38 


< 0.62 


< 1.67 


55 


SR21 


< 19.59 


< 24.43 


< 87.39 


56 


RX.J1852.3-3700 


< 1.26 


0.80 ± 0.11 


< 0.89 



Table 3. [Ne II] /[Ne III] Line Flux Ratio Literature Data 



Literature data ID Name [Ne II] /[Ne III] flux ratio Ref. 



1 


IRAS 03446+3254 


< 2.77 


1 


2 


IRAS 08267-3336 


< 3.5 


1 


3 


Ced 110 IRS 6 


< 2.38 


1 


4 


VW Cha 


< 5.09 


1 


5 


XX Cha 


< 3.89 


1 


6 


T Cha 


< 16.00 


1 


7 


Sz 73 


< 4.57 


1 


8 


IM Lup 


< 2.26 


1 


9 


V853 Oph 


< 4.58 


1 


10 


IRS 60 


< 2.15 


1 


11 


Haro 117 


< 2.67 


1 


12 


SSTc2d J182928.2+02257 


< 5.33 


1 


13 


EC 74 


< 3.13 


1 


14 


EC 92 


< 2.63 


1 


15 


Sz 102 


15.65 ± 4.8 


1 


16 


RX Jllll. 7-7620 


< 1.7 


2 


17 


[PZ99] J161411.0-230536 


< 1.32 


2 


18 


RX J1842.9-3532 


< 1.19 


2 


19 


RX J1852. 3-3700 


< 4.80 


2 


20 


TW Hya 


22.24 ± 16.9 


3 


21 


WL5/GY246 


3.75 ± 5.3 


4 


22 


DoAr25/GY17 


< 2.04 


4 


23 


WL12/GY111 


< 2.12 


4 


24 


WL10/GY211 


< 1.95 


4 


25 


WL20/GY240 


< 4.07 


4 


26 


IRS37/GY244 


< 4.34 


4 


27 


IRS43/GY265 


< 7.28 


4 


28 


IRS44/GY269 


< 1.04 


4 


29 


IRS45/GY273 


< 3.54 


4 


30 


IRS47/GY279 


< 1.80 


4 


31 


DM Tau 


< 0.67 


5 


32 


Coku Tau-4 


< 4.0 


5 



-44- 



Table 3 — Continued 



Literature data ID Name [Ne II] /[Nc III] flux ratio Rcf. 



33 GM Aur < 1.5 5 



Najit^_e^al 



Baldovin-Saavedra et al 



Not e. — References: (l ) lLahuis et al.N2007l'l : (2 
I 2007l'l : (3) 



J gold) : (4) 



Pascucci et al. 



Flaccomio et al 



1 20091 '): (5) 



i 201lh 



